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It has recently been shown that a non-Hermitian Hamihonian H possessing an 
unbroken VT symmetry (i) has a real spectrum that is bounded below, and (ii) 
", defines a unitary theory of quantum mechanics with positive norm. The proof of 

I unitarity requires a linear operator C, which was originally defined as a sum over the 

eigenf unctions of H. However, using this definition to calculate C is cumbersome in 
quantum mechanics and impossible in quantum field theory. An alternative method 



o 

(N 



D . is devised here for calculating C directly in terms of the operator dynamical variables 



m 
(N 



of the quantum theory. This new method is general and applies to a variety of quan- 
tum mechanical systems having several degrees of freedom. More importantly, this 
method is used to calculate the C operator in quantum field theory. The C operator 
is a new time-independent observable in PT-symmetric quantum field theory. 
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2 ■ I- INTRODUCTION 

o 

In 1998 numerical and perturbative methods were used to establish the surprising result 
that the class of non-Hermitian Hamiltonians 

H = p^ + x\ixy (e>0) (1) 

^ ■ 

has a positive real spectrum In Ref. |l| it was conjectured that the spectral positivity was 
I associated with the space-time reflection symmetry (VT symmetry) of the Hamiltonian. The 
■ Hamiltonian if in is VT symmetric because under parity reflection V we have x —x 
^ and p — >■ —p and under time reversal T we have x — x, p — > — p, and i — > —i. Other 

PT-symmetric quantum mechanical models have been examined 0, 0, 0, S 0] , and a proof 
of the positivity of the spectrum of if in ((T)) was subsequently given by Dorey et al. 'T]. 

The discovery that the spectra of many PT-symmetric Hamiltonians are real and positive 
raised a fundamental question: Does a non-Hermitian Hamiltonian such as H in ^ define 
a consistent unitary theory of quantum mechanics, or is the positivity of the spectrum 
merely an intriguing mathematical property of special classes of complex Sturm-Liouville 
eigenvalue problems? To answer this question it is necessary to know whether the Hilbert 
space on which the Hamiltonian acts has an inner product associated with a positive norm. 
Furthermore, it is necessary to determine whether the dynamical time evolution induced by 
such a Hamiltonian is unitary; that is, whether the norm is preserved in time. 

Recently, a definitive answer to this question was found [1, |^. For a complex non- 
Hermitian Hamiltonian having an unbroken VT symmetry, a linear operator C that com- 
mutes with both H and VT can be constructed. The inner product with respect to CVT 
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conjugation satisfies the requirements for tlie tlieory defined by H to have a Hilbert space 
with a positive norm and to be a consistent unitary theory of quantum mechanics. [The 
term unbroken VT symmetry means that every eigenfunction of H is also an eigenfunction 
of the VT operator. This condition guarantees that the eigenvalues of H are real. The 
Hamiltonian in has an unbroken VT symmetry for all real e > 0.] 

We emphasize that in a conventional quantum theory the inner product is formulated 
with respect to ordinary Dirac Hermitian conjugation (complex conjugate and transpose). 
Unlike conventional quantum theory, the inner product for a quantum theory defined by a 
non-Hermitian "PT-symmetric Hamiltonian depends on the Hamiltonian itself and is thus 
determined dynamically. One can view this new kind of quantum theory as a "bootstrap" 
theory because one must solve for the eigenstates of H before knowing what the Hilbert 
space and the associated inner product of the theory are. The Hilbert space and inner 
product are then determined by these eigenstates. 

The key breakthrough in understanding these novel non-Hermitian quantum theories was 
the discovery of the operator C j^. This operator possesses three crucial properties. First, 
it commutes with the space-time reflection operator VT, 

[C,PT]=0, (2) 

although C does not commute with V or T separately. Second, the square of C is the identity, 

C = 1, (3) 

which allows us to interpret C as a reflection operator. Third, C commutes with H, 

[C,H] = 0, (4) 

and thus is time independent. To summarize, C is a new time-independent PT-symmetric 
reflection operator. 

The question now is how to construct C for a given H. In Refs. it was shown how 

to express the C operator in coordinate space as a sum over the appropriately normalized 
eigenfunctions of the Hamiltonian H. These eigenfunctions satisfy 

H(t)n{x) = EnCpnix), (5) 

and, without loss of generality, their overall phases are chosen so that 

VT^x) =(j)n{x). (6) 
With this choice of phase, the eigenfunctions are then normalized according to 

dx [M^)f = (-1)". (7) 

c 

The contour of integration C is described in detail in Ref. For the quantum mechanical 
theories discussed in this paper, all of which have a cubic interaction term, the contour C 
can be taken to lie along the real-a; axis. 

In terms of the eigenfunctions defined above, the statement of completeness for a theory 
described by a non-Hermitian PT-symmetric Hamiltonian reads 0] 

J2i-^rM^)My) = Six - y) (8) 
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for real x and y. The coordinate-space representation of C is 

C{x, y) = Y^ (j)n{x)(f)n{y). (9) 

n 

Only a non-Hermitian PT-symmetric Hamiltonian possesses a C operator distinct from 
the parity operator V. Indeed, if one evaluates the summation Q for a "PT-symmetric 
Hamiltonian that is also Hermitian, the result is V, which in coordinate space is 6{x + y). 
The coordinate-space formalism using Q has been applied successfully to 

H = ^p^ + i/i^x^ + iex\ (10) 

and C was constructed perturbatively to order e'^ pjj . This formalism has also been applied 
to calculate C to order e for the complex Henon-Heiles Hamiltonian 

h = 1{pI+pI) + 1 + y') + '^^'y^ (11) 

which has two degrees of freedom, and for the Hamiltonian 

H = 1 {px^ +pI+pI) + 1 {x^ + y^ + z^) + iexyz, (12) 

which has three degrees of freedom jl3| . 

Calculating the operator C by direct evaluation of the sum in is difficult in quantum 
mechanics because it is necessary to determine all the eigenf unctions of H . Such a procedure 
cannot be used at all in quantum field theory because there is no simple analog of the 
Schrodinger eigenvalue problem © and its associated coordinate-space eigenfunctions. 

In this paper we devise an elementary operator technique for calculating C for the impor- 
tant case of quantum theories having cubic interactions, and we demonstrate that our new 
method readily generalizes from quantum mechanics to quantum field theory. In Sec. |H] we 
introduce a general operator representation for C of the form e'^^^'^^V, where x and p are 
the dynamical variables. This representation is especially convenient for incorporating the 
three requirements © ^ ©• In Sec. II I II we calculate C to seventh order in powers of e for 
the Hamiltonian pO|) using this operator technique. Then, in Sec. II VI we calculate C for the 
Hamiltonians (fTT|) and (fT^ to order e^. In SecElwe apply operator methods to calculate C 
for the massless Hamiltonian H = + ix^. We derive recursion relations for the operator 
representation of C in Sec. IVII In Sec. I VIII we calculate C to order for the self- interacting 
scalar quantum field theory described by the Hamiltonian 

H = y"c/^x|i7r2(x,t) + i[Vx^(x,t)]2 + i/iV'(x,t)+ie(/?3(x,t)} (13) 

in [D + l)-dimensional Minkowski spacetime. In Sec. IVIIII we calculate C for cubic scalar 
quantum field theories with interactions of the form ie(p\ip2 and ie(fi(p2f3- An alternative 
perturbative calculation of C for an ie^p^ quantum field theory using diagrammatic and 
combinatoric methods is given in Appendix ^ Some concluding remarks are in Sec. IIXI 

The principal accomplishment of this paper is the derivation in Sees. IVIII and IVIIII and 
Appendix ^of the C operator for cubic quantum field theories. Cubic quantum field theories, 
such as that in (jl3|) . are not just of mathematical interest. Such theories emerge in the study 
of Reggeon field theory and in the analysis of the Lee- Yang edge singularity For 
these quantum field theories the operator C is a new conserved quantity. Knowing how 
to calculate this operator is crucial because it is necessary to have C in order to construct 
observables and to evaluate matrix elements of field operators. Our calculation of C is a 
major step in our ongoing program to obtain new physical models by extending conventional 
quantum mechanics and quantum field theory into the complex domain. 
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II. GENERAL FORM FOR THE OPERATOR C 

To prepare for calculating C we show in this section that it is advantageous to represent 
C as a product of the exponential of a Hermitian operator Q and the parity operator V: 

C = e^^^'^'^P. (14) 



This representation was first noticed in Ref. [11 



A. Previous Work on Calculating C 

The objective of the investigation in Ref. [ll[ was to use perturbative methods to calculate 
C for the Hamiltonian H = ^p^ + ^x"^ + iex^, where e is treated as a small parameter. In 
Ref. [llj the operator C was obtained to third order in e in coordinate space. The procedure 
was as follows: First, the Schrodinger equation 

- |0"(a;) + |a;^0„(x) + iex^0„(x) = E„0„(x) (15) 

was solved for the energies En and for the wave functions (f)n{x) as Rayleigh- Schrodinger 
perturbation series in powers of e. The series for (/)„(x) has the form 

0„(x) = r \ rr ^"'"'' [Hn{x) - tAn{x)e - Bn{x)^ + za,(x)e^ + ■■■], (16) 

where Hn{x) is the nth Hermite polynomial and A„(x), Bn{x), and C„(x) are polynomials 
in X of degree n + 3, n + 6, and n + 9, respectively. These polynomials were expressed as 
linear combinations of Hermite polynomials. The value of a„, 

a„ = 1 + ^(2n + l)(82n2 + 82n + 87)6^ + O(e^), (17) 

ensures that the eigenfunctions are normalized according to (ix0^(x) = (—1)", as in ((7j). 
The factor i"' in (fT^ is included to satisfy the requirement in (jH)) that VT(j)n{x) = 0„(x). 

Finally, (fT6|) was substituted into (jnj and the summation over n was performed to obtain 
the operator C{x,y) to order e^: 

C{x, y) = 1 1 - e - 2xyp) + [§/ - fxy/ + {2xY ~ 12)p'] 

-(8x^7/2 - 64)p] + 0(e^)|5(a; + y), (18) 

where p = —id/dx. This expression is complicated, but it was observed that it simplifies 
considerably when the expression in curly brackets is rewritten in exponential form: 

C{x, y) = e'^'+''^'+-6{x + y) + O(e^). (19) 

In this form the differential operators Qi and Q3 are simply 

Qi = - 2xpx, 

Q3 = ^fxp'x + Sx^px'^ -Up. (20) 

The main features of the exponential representation (fT^ are that only odd powers of e 
appear in the exponent, the coefficients are all real, and the derivative operators act on the 
parity operator S{x + y). Also, e^*^^"^^ is Hermitian. 
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B. New Approach to Calculating C 

The perturbative calculation described above suggests that a simpler and more direct way 
to calculate C is to seek an operator representation of it in the form e'^^^'^^P, where Q{x,p) 
is a Hermitian function of the operators x and p. We will show that Q{x,p) can be found by 
solving elementary operator equations and that it is not necessary to find the eigenfunctions 
to determine Q. Thus, the technique introduced in this paper immediately generalizes to 
quantum field theory. To find the operator equations satisfied by Q we substitute C = e^V 
into the three equations (0) - (jH). The details are described in Sec. IIIII 

We claim that the representation C = e^V is general. Let us illustrate this simple 
representation for C in two elementary cases: First, consider the shifted harmonic oscillator 

H = \p' + \x^ + iex. (21) 

This Hamiltonian has an unbroken VT symmetry for all real e. Its eigenvalues En = 
+ I + are all real. The C operator for this theory is given exactly by C = e'^V, where 
Q = —ep. Note that in the limit e ^ 0, where the Hamiltonian becomes Hermitian, C 
becomes identical with V. 

As a second example, consider the non-Hermitian 2x2 matrix Hamiltonian 

ff= ('■:%;«). (22) 

which was discussed in Ref. jlo|. This Hamiltonian is VT symmetric, where V is the Pauli 
matrix ai = (H) and T is complex conjugation. This Hamiltonian has an unbroken VT 
symmetry when > sin^ 9. The C operator in the unbroken region is 



cos a 



i sin a 1 
1 —i sin a 



(23) 



where sin a = (r/s) sin 6. Our new way to express C is to rewrite it in the form C = e^V. 
Thus, the Hermitian operator Q has the form 

Q = la2lni-—. , 24 

^ \ 1 + sm a / 

where (T2 = (^ q*)- Again, observe that in the limit 6^0, where the Hamiltonian becomes 
Hermitian, the C operator becomes identical with V. 

Note that V, which is given in coordinate space as 6{x + y), can be expressed in terms of 
the fundamental operators x and p as 

V{x,p) = exp [|OT(p2 + - 1)] . (25) 

To show that the parity operator satisfies VxV^^ = —x and VpV^^ = —p, we define the 
operator- valued functions /(r) and g^r) as 

/(r) = e^-(p'+^')xe-^"(p'+^'), g^r) = eMP^+^')pe-^"(p'+^'). (26) 

Differentiating /(r) and g^r) once gives gives /'(r) = 2g{T) and g'ir) = — 2/(r). A second 
differentiation then leads to the differential equations /"(t) = 4/(r) and g"{T) = 4(yf(r). 
The solutions to these equations satisfying the initial conditions /(O) = x and g{0) = p are 

/(r) = a; cos(2r) — zpsin(2r), 5'(r) = pcos(2r) — sin(2r). (27) 
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Setting T = iyr, we get /(r) = —x and g{T) = —p, which estabhshes that the operator V 
defined in fj25p indeed has the properties of a parity refiection operator. Specifically, P is a 
unitary operator that generates a rotation by vr in the {x,p) plane. Another application of 
V gives a rotation by 27r in the {x,p) plane. Hence = 1. This procedure determines V 
up to an additive phase. It is conventional to choose the phase to be — |vr, as in ()25|1 . 

III. CUBIC OSCILLATOR WITH ONE DEGREE OF FREEDOM 

Having shown that (|n|l is a natural way to represent the operator C, we now demonstrate 
how to use this ansatz to calculate C for the Hamiltonian The procedure is to impose 
the three conditions Q - Q in turn on C = e^'^^'^^V and thereby to determine the operator- 
valued function Q{x,p). 

First, we substitute ()14|) into the condition Q to obtain 

gQ(x,p) _ 'pq-QQix,p)'pq- _ ^Q{-x,p)^ 

from which we conclude that Q{x,p) is an even function of x. Second, we substitute (fT^ 
into the condition (jS)) and find that 

gQ(x,p)-pgQ(x,p)p _ ^Q{x,p)^Q(-x-p) _ -^^ 

which implies that Q{x,p) = —Q{—x, —p). Since we already know that Q{x,p) is an even 
function of x, we conclude that it is also an odd function of p. 

The remaining condition Q to be imposed is that the operator C commutes with H. 
Substituting C = e^^^'^^P into ©, we get e'3(^'P)[P, i/] + [e^^'^'P\H]V = 0. All of the 
Hamiltonians H considered in this paper can be expressed in the form H = Ho + eHi, 
where Ho is a free field theory (harmonic oscillator) Hamiltonian that commutes with the 
parity operator V, and Hi represents the interaction. For example, for the Hamiltonian 
H = + + iex^, Hq = + and Hi = ix^. Then, this condition reads 

e'3(^'P)[P, Hi] + [e^(^'P\ H]V = 0. (28) 

Next, we observe that if the interaction is cubic, then Hi is oc?c? under parity refiection; that is. 
Hi anticommutes with V. Hence, for quantum theories with cubic interaction Hamiltonians, 
(OHj) reduces to 

eeQ(^',P)^^ = [eQ{^,P)^^]_ (29) 

We note that the structure in (fT^ is quite general; for all cubic Hamiltonians, Q{x,p) 
may be expanded as a series in odd powers of e: 

Q{x,p) = eQi{x,p) + e^Qsix^p) + e^Q^ix^p) + ... . (30) 

In quantum field theory we will interpret the series coefficients Q2n+i as interaction vertices 
(form factors) of 2n + 3 powers of the quantum fields. 

Substituting the expansion in (jHUj) into the exponential e'^^^'P\ we get 



^Qix,p) ^ j^^^^p^ ^ ^ ^ Ri{x,p)e + R2{x,p)e^ + R3{x,p)e'^ + Ri{x,py + . . . , (31) 
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where 

Ri = Qi, 

R2 = iQi, 

R3 = Qs + 

-R4 = ^{QijQs} + j^Qi^ 

R5 = Q5 + l{{Qi,Q3} + Q1Q3Q1) + j^Qi, 

Re = UQI + {^1' Q^}) + UiQl Q^} + {^1' Q1Q3Q1}) + WoQl^ (32) 

and so on. Here, {X, Y} = XY + YX denotes the anticommutator. 

We now substitute (jH^ into (jSHl), collect the coefficients of like powers of for n = 
1, 2, 3, . . ., and obtain a sequence of equations of the general form 

e^: [Ho,Rn] = -{H^,Rr,^^} (n > 1), (33) 

where Rq = 1. 

The equations in (j33j) can be solved systematically for the operator-valued functions 
Qn{x,p) {n = 1, 3, 5, . . .) subject to the symmetry constraints that ensure the conditions ^ 
and (jni). Note that the coefficients of even powers of e contain no additional information 
because the equation arising from the coefficient of e^" can be derived from the equations 
arising from the coefficients of e^"~^, e^""^, . . . , e. This observation leads to a more effective 
way to express the conditions in (jHSJ- The ffist four equations read 

[Ho,Qi] = -2Hi, 

[Ho,Q3] = -llQulQuH,]], 

[Ho,Q6] = ^,[Qu[Qi,[QiAQuM]]-li[QiAQ3,H,]] + [Q,,[Q,,H,]]), 

[Ho, Q7] = T5T20 [^1' [^1' [Qi' [Qi' [^1' [^1' ^1]]]]]] 

-sio ([^1' [^1' [^1' [^3, H^]]]] + [Qi, [Q^[Q,, [Q,, H,]]]] 

+ [Qu [Q3, [Qi, [QuH,]]]] + [Q3, [QiiQu [QuH,]]]]) 

+1 i[Q^, [Q,,H,]] + [Q,, [Q,,H,]]) l[Q,, [Q^, H,]]. (34) 

We now show how to solve these equations for the Hamiltonian in (fTUj) . for which Hq = 
+ and Hi = ix^. The procedure is to substitute the most general polynomial form 

for Qn using arbitrary coefficients and then to solve for these coefficients. For example, to 
solve the ffist of the equations in (jH^ . [Ho,Qi] = —2ix^, we take as an ansatz for Qi the 
most general Hermitian cubic polynomial that is even in x and odd in p: 

Qi{x,p) = Mp^ + Nxpx, (35) 

where M and are undetermined coefficients. The operator equation for Qi is satisfied if 

M = and = -2fx-^. (36) 

It is straightforward, though somewhat tedious, to continue this process. In order to 
present the solutions for Qnix, p) {n> 1), it is convenient to introduce the following notation: 
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Let Sm,n represent the totally symmetrized sum over all terms containing m factors of p and 
n factors of x. For example, 



'S'0,0 = 1; 

So,3 = x"^, 

51.1 = l{xp + px) , 

'S'1,2 = I {x^P + xpx + px^) , 

Sz,i = \ {xp^ + pxp^ + p^xp + p^x) , 

82.2 = \ {p^x^ + x^p^ + pxpx + xpxp + px^p + xp^x) . (37) 



The properties of the operators Sm,n are summarized in Ref. [17[. One useful property is 
that Sm,n can be expressed in Weyl-ordered form in two ways: 

fc=0 fc=0 

We have solved the equations in and have found Qi, Q3, Q^, and Qy in closed form. 
In terms of the symmetrized operators Sm,n the functions Qn are 

Qi = -|/^"V - 2/i-'5i,2, 

Qs = ^fi-'V + ffi-'Ss,2 + 8/i-^^i,4 - I2fi-'p, 

n — 320,-16„7 544, -14 Q 512,-12 r- 
V5 — 3-/^ P 3-/^ '-'5,2 3-/^ '-'3,4 

fiA/z-lO'? I 24 736 , ,-14 3 , 6368,, -12 r. 
— 04/i Di_6 H 45—/^ P H 15-^^ '-'1,2, 

n _ 553 984 ,,-22^9 I 97 792 ,, -20 c 1 377344 ,, -18 c 
V7 - 315 P + -35"^ "^7,2 + 105 '^5,4 

I 721024 ,,-16 c I 1792 ,, -14c 2 209 024 ,,-20^5 
+ 315 '-'3,6 H 3— A* '-'1,8 105 — ^ ^ 

2875648 ,,-18 c 390 336 ,,-16 c 1 46976 , ,-18^ /qQ\ 
A* '-'3,2 35 — A* '-'1,4 H 5 — P- \^^) 



105 ^ ^■^'^ 35 

Combining (HH), (jHUj) . and (jH^ . we obtain an explicit perturbative expansion of C in terms 
of the fundamental operators x and p, correct to order e^. 

To summarize, using the ansatz (|T1| we are able to calculate the C operator to very high 
order in perturbation theory. We are able to perform this calculation because this ansatz 
obviates the necessity of calculating the wave functions (t)n{.x). The calculation bears a strong 
resemblance to WKB theory. The ansatz used in performing a semiclassical calculation is 
also an exponential of a power series. The advantage of using WKB to calculate the energy 
eigenvalues is that to all orders in powers of h it is possible to construct a system of equations 
like those in ()34|) that determine the energies, and it is never necessary to calculate the wave 
function Furthermore, only the even terms in the WKB series are needed to determine 
the energy eigenvalues. The odd terms in the series drop out of the calculation and provide no 
information about the eigenvalues [3]. The difference between a conventional WKB series 
and the series representation for Q is that the first term in a WKB series is proportional to 
h^^ while the series expansion for Q{x,p) contains only positive powers of e. Based on the 
results in Ref. however, we believe that for a PT-symmetric —ex^ theory, the first term 
in the expansion of Q{x,p) is proportional to e~^. We plan to discuss quartic PT-symmetric 
theories in a future paper. 
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IV. CUBIC OSCILLATORS WITH SEVERAL DEGREES OF FREEDOM 

In this section we extend the operator techniques used in Sec. IIIII to systems having 
two and three dynamical degrees of freedom. Specifically, we generalize the perturbative 
procedure for calculating the C operator for the Hamiltonian in and use it to calculate 
C for the Hamiltonians in (jllj) and ()12j) . 

Let us first consider H in which has two degrees of freedom. We write this Hamil- 
tonian in the form H = Hq + eHi, where Hq = + + + \y'^ and Hi = ix^y. For 
these operators we need to solve the system of equations in (jMj) for the unknown operators 
Qii Q3, and so on. To simplify the calculation we generalize slightly the notation in (jHTj) for 
totally symmetric operators; to wit, we continue to use Sm,n to represent a totally symmetric 
product of m factors of p and n factors of x but we use Tm,n to represent a totally symmetric 
product of m factors of q and n factors of y. For example, 

^1,1 = U^y + yq), 

Ti,2 = I {y'^q + yqy + qy'^) = yqy, 

Ts,! = I {yq^ + qyq^ + q'^yq + q^y) , 

7^2,1 = \ [q^y + qyq + yq^) = qyq- (40) 

To solve [Hq, Qi] = —2Hi, the first equation in (|34p . we seek a Hermitian cubic polynomial 
in the variables x, y, p, and q. This polynomial must be even in the coordinate variables 
and odd in the momentum variables, and it must be able to yield Hi when commuted with 
Hq. We therefore introduce the ansatz 

Qi{x,y,p,q) = Mp^q + NiSi,iy + N2X^q. (41) 

We substitute this ansatz into the commutator and determine the unknown constants M, 
Ni, and N2 by solving three linear equations. The result is 

M = -|, iVi = -|, N2 = -l (42) 

Next, we turn to the second of the equations in (jSlI), [-f^o^Qs] = ^^IQi, [Qi,Hi]], and 
evaluate its right side. The resulting equation for then reads 

[^0, Qs] = ^§f (4x^2/ + AxY + 8p'T2,i - 4x^^,1 + ASi,^q 

+8^i,iri,2 + SSs,iq - 3y) . (43) 

We now must construct the most general Hermitian fifth-degree polynomial in the variables 
X, y, p, and q that is even in the coordinate variables, odd in the momentum variables, and 
has the terms needed to produce the right side of this commutation relation: 

Q3{x,y,p,q) = aip^q^ + a2P^q + a^Si^iT2,i + aip^Ti^2 + a'bSz,iy + a'%S2,2q 

+a-jx'^q^ + asx'^Ti^2 + agSi^iy^ + aioSi^^y + anx^q + auq- (44) 

Substituting this ansatz into (jlHj) . we obtain twelve simultaneous linear equations for the 
unknown coefficients a^, whose solution is 
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This completes the calculation of the operator C to third order in e for H in 

The attractive feature of the calculational procedure described here is that it is utterly 
routine and works in every order of perturbation theory. In contrast, the technique used 
in Ref. [l^ to calculate C becomes hopelessly difficult beyond first order in powers of e 
because the technique used earher requires that one calculate all of the energy eigenstates 
(f>n{x, y) perturbatively for all n. These eigenstates must then be substituted directly into the 
summation in that defines the operator C. This calculation is difficult because beyond 
leading order in e one encounters the challenging problems associated with degenerate energy 
levels. (There are no degenerate energy levels for Hamiltonians having just one degree of 
freedom.) Of course, for any given n, there is a well-defined procedure for calculating the 
eigenstate to any order in powers of e. However, this procedure depends on the value of n and, 
as a result, the calculation becomes extremely complicated. The new method of calculation 
presented here works because it is no longer necessary to calculate the eigenfunctions. Thus, 
the difficulties associated with degeneracy are circumvented. 

We turn next to the case of a cubic oscillator having three degrees of freedom. We express 
the Hamiltonian in (fT^ in the form H = H^+eHi, where Hq = +\q'^ ^-^r"^ +\x'^ +\y'^ +^ ""^ 



z 



represents a harmonic oscillator Hamiltonian having three degrees of freedom and Hi = ixyz 
is a non-Hermitian PT-symmetric interaction term. 

To solve [Hq, Qi] = —2Hi, the first of the equations in (jH^ . we must construct the most 
general Hermitian cubic polynomial in the variables x, y, z, p, q, and r that is even in the 
coordinate variables, odd in the momentum variables, and has the terms needed to yield Hi 
on the right side of this commutation relation: 

Qi{x,y, z,p,q,r) = Mpqr + N{yzp + xzq + xyr). (46) 

We then substitute this ansatz into the commutator and determine the unknown constants 
M and by solving two linear equations. The result is 

M = -| and = -f . 

To solve the second of the equations in (jH^ . [Hq, Qs] = —\[Qii [<5i, H-\\], we evaluate its 
right side. The resulting equation for Q3 then reads 

[Hq, Qs] = xyz {x^ + y"^ + z^) + {pz + rx)yqy + {py + qx)zrz + {qz + ry)xpx 

— {pxpyz + qyqxz + rzrxy) + 2{rqpxp + rpqyq + pqrzr) . (47) 

The most general Hermitian fifth-degree polynomial in the variables x, y, z, p, q, and r 
that is even in the coordinate variables, odd in the momentum variables, and has the terms 
needed to produce the right side of this commutation relation is: 

Qsix, y, z,p, q, r) = di {p^qr + q^pr + r^gp) + d2[pxp{yr + zq) + qyq{xr + zp) 

+rzr{xq + yp)] + d^^xpxqr + ypypr + zrzpq) + d^^xpxyz + yqyxz 

+zrzxy) + d^ [x'^(?/r + zq) + y^{xr + zp) + z^{xq + yp)] 

+dQ {p^yz + q^xz + r^xy) . (48) 

Substituting this ansatz into (jlHjl . we obtain six simultaneous equations for the unknown 
coefficients dn-, whose solution is 
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This completes the calculation of the operator C to third order in e for H in ()12|) . This 
calculation is simpler than that for the Hamiltonian in because there is symmetry 
under the interchange of pairs of dynamical variables, such as {x,p) ^ {y,(l)- 

V. REPRESENTATION OF C FOR THE MASSLESS CASE 

In this section we examine the massless (strong-coupling) limit of the operator C for H 
in (|Tn|l . The massless limit /i — of the massive theory is especially interesting because 
this limiting case is singular. We will see that as the mass parameter /i tends to zero, the 
perturbation series representation for Q in C = e'^P ceases to exist and an entirely new 
nonpolynomial representation for Q emerges. 

Negative powers of /i in fl39|) are required for dimensional consistency. As a result, each 
of the terms in these perturbation series coefficients becomes singular in the massless limit 
/i — > 0. (Note that the dimensionless perturbation expansion parameter is e/i~^/^, and thus 
the massless limit of the theory is equivalent to the strong-coupling limit e — > oo.) 

To find the C operator for the massless theory it is necessary to return to the sequence 
of operator equations in (j34p and to look for new solutions for the special case in which 
Ho = The first of these equations reads 

=-2tx\ (50) 

However, an examination of (jKn|) reveals that it is no longer possible to find a solution in the 
form of a polynomial in the operators p and x. The situation here is quite similar to that 
considered in Ref. in which the objective was to calculate the time operator in quantum 
mechanics. In the case of the time operator it was shown that one must generalize the 
symmetric operators Sm,n from the positive integers to the negative integers. Specifically, if 
m is nonnegative, then n may be negative, and if n is nonnegative, then m may be negative. 
(It is not possible for both m and n to be negative.) We can display these generalized 
symmetric operators in Weyl-ordered form [see p8|l ]. For example, 

l(x'-+'-x], 

2 \ p p J' 
J_ 

l[^x' + 2x^x + x'^), 

l{x^^ + Sx'jjx + Sxjjx^ + ^x^) . (51) 
An exact, dimensionally consistent operator solution to ()50|) is 

Qi = i5_i,4 + a5„5,o = 3^ {x^l + ^x'lx + Qx^y + Ax^ + ^a;^) + a^, (52) 

where a is an arbitrary number. This solution has the required symmetry properties; to 
wit, it is odd in p and even in x. Also, it has the same dimensions as Qi in ()39|) . 

The solution to the second operator equation in [|p^,Q3] = — [<5i, '^a^'^]], is 

— ^'S'_5,io — ^S^-j^s + + 20a) S'_9,6 + (m ^'^) '^'-11,4 

+ (- W - ^« + ^-13,2 + /3^-i5,o, (53) 



'S'-3,0 
S-2,2 

S-2,3 
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where (3 is a new arbitrary constant. Again, observe that this solution exhibits the required 
symmetry properties. 

A notable feature of the solutions for Qi and Q3 in and is that they are not 
unique. Each of these solutions contains an arbitrary constant multiplying a negative odd- 
integer power of p. There is no obvious way to determine the values of the constants a and 
j3. These terms arise because in the massless case the Hamiltonian Hq is a function oip only. 
In general, one can add an arbitrary multiple of p-iO"-^ to the solution for Q2n+i because it 
is odd in p and is dimensionally consistent. In the massive case, where Hq = + 
there is no such ambiguity because adding an arbitrary function of Hq to Qn would violate 
the symmetry requirement that Qn be odd in p. 



VI. PRODUCT REPRESENTATION OF C AND DERIVATION OF 

RECURSION RELATIONS 

In this section we investigate the product representation of the operator C that was defined 
in ()31|): namely C{x,p) = R{x,p)V. At first, it may not seem worthwhile to reconsider the 
product representation because it lacks the advantages of the exponential representation C = 
e^V introduced in (|14j). Recall that we argued in Sec. Othat the exponential representation 
is convenient because it incorporates the requirements (0) and Q as elementary symmetry 
conditions on Q{x,p): Q{x,p) = Q{—x,p) and Q{x,p) = —Q{x,—p). Furthermore, we 
showed that the exponential representation of C in and piU) is much simpler than the 
product representation C{x,p) = R{x,p)V in (fTH|) . However, as we demonstrate here, the 
product representation has the advantage that it can be used to construct a recursive formula 
for the perturbation coefficients. 

The function R{x,p) in the product representation C{x,p) = R{x,p)V incorporates the 
requirement in Q as 

R{x,p) = R{-x,p). (54) 

Thus, R{x,p) is an even function of x. However, the requirement in © translates into a 
complicated nonlinear condition on R: 

R{x,p)R{x,-p) = 1. (55) 

We will return to this condition later. 

The advantage of the product representation is that it translates the requirement in Q 
into a linear difference equation. To obtain this difference equation we use the operators 
Sm,n in (IHTjl . (Recall that Sm,n is a totally symmetric combination of products of m factors 
of p and n factors of x.) It was shown in Ref. jl^ that the operators Sm,n are complete in 
the sense that any operator may be represented as a linear combination of these symmetric 
operators. This allows us to represent R{x,p) as the infinite linear combination 

C^m,nSm,ni (^^) 

m n 

where am,n are numerical coefficients to be determined. Substituting C{x,p) = R{x,p)V 
into (0} then gives the condition 



[R,Ho]=e{R,Hi} 



(57) 
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We now substitute R in (I56p into the condition ()57p and use the commutation and anti- 
commutation relations lO] 

{Sni,n,X^] = -lm{m- l)Sm~l,n+l- (58) 

For the Hamiltonian in ()10|) we obtain the hnear recursion relation 

n am-2,n - fi'^m a,rn,n-2 = e [-|m(m + 1) a„+i^„_2 + 2ara-l,n-4\ ■ (59) 

The boundary conditions on this partial difference equation must be chosen so that the 
nonlinear constraint (j^^ is satisfied. In the massive case ao,o = 1 and am,n vanishes if either 
m < or n < 0. In the massless case we again have ao,o = 1, but now am,n vanishes if either 
m > or < 0. 

One approach to solving this equation is to introduce a generating function g{s,t): 

^(s,t) = ^^«^,„s™t". (60) 

m n 

For the massive case [the Hamiltonian in (fTUj) with /i 7^ 0] the summation is taken over 
nonnegative values of m and n. However, for the massless case (/i = 0) the summation must 
be taken over nonpositive values of m and nonnegative values of n. 

We then multiply ()59p by s"^~^t"^~^ and rewrite the result in the form 

^'Ol V^m-2,nS t ) - jj, t— (am,„-2S t ) 



rj2 

3+ (n, „m+l.n-2\ , 9.3 m-1 .n-4 

'2^'g^ ['^rn+l,n-2S t ) + Zt am-l,n-AS t 



(61) 



Summing over m and n and using ()60|). we obtain the partial differential equation 

sgt - fiHgs = e (- ft + 2t^g) , (62) 

where subscripts indicate partial differentiation. A Fourier transform from the s variable 
to the r variable converts this differential equation into the Goursat problem (recall that a 
Goursat problem involves a wave equation written in light-cone variables) 

9rt = {fJ''^rt + ^ier^t + 2iet^) g. (63) 

This is an extremely interesting equation that merits further study. We plan to give a 
detailed analysis of this partial differential equation and of the partial difference equation 
dnni) in a future paper. 



VII. SCALAR QUANTUM FIELD THEORY WITH CUBIC 

SELF-INTERACTION 

This section extends the operator techniques introduced in Sec. IIIII to quantum field 
theory. Consider the quantum field theory described by the Hamiltonian (fT!?|l in {D + 1)- 
dimensional Minkowski space-time. This Hamiltonian has the form H = Hq + eHi, where 

Ho= /d^a;{i7r2(x,t) + i[Vx<^(x,t)]2 + i^V'(x,t)}, = i [d^'x ip%x,t). (64) 
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The integrals above are performed in the spatial variable x, which lies in R^. In the 
following we use J (ix = J d^x to represent the integration in MP . The field variables satisfy 
the equal-time canonical commutation relation [<^(x, t), 7r(y, t)] = i(5(x — y). 

The parity operator is given formally by "P = exp (iivr / dx [¥?2(x,t) + 7r2(x,t) - 1]). As 
in quantum mechanics, where the operators x and p change sign under parity reflection, we 
assume that the fields are pseudoscalars and that they also change sign under V: 

P(^(x, t)V = -<^(-x, t), P7r(x, t)V = -7r(-x, t). (65) 

Following the approach in Sec. IIIIl we express C in the form C = e^Qi+'^^Qs+-'p ^ where 
Q2n+i {n = 0, 1, 2, . . .) are real functionals of the field variables (/^(x, t) and 7r(x, t). To find 
Qi it is necessary to solve the first of the operator equations in (jH^ : 

y"rfx(i7r2(x,t) + i/iV'(x,t)-i^(x,t)V,V(x,t)),Qi =-2t Jd^ip\^,t), (66) 

where we have integrated by parts: J (ix(VxV5)^ = — / d^KipV^ip. We define the inverse 
Green's function G"^ by G^^ = {fi^ - V^)6{^ - y), so that G^y = (/i^ - V2)-i(5(x - y) and 
fdzG~^Gzy = 5(x — y). Thus, the commutator condition (jUHjl reads 

i y"c/x7r2(x,t) + i J dxdyipix,t)G^^ipiy,t),Qi = ~2i jdxip%x,t). (67) 

Equation ()67|1 states that when the operator Qi is commuted with quadratic structures of 
the form 7r^(x, t) and ip{:s.,t)ip{y,t), it must produce the cubic term (y9^(x, t). Furthermore, 
the symmetry requirements on Qi that arise from ^ and (jH)) imply that Qi is an even 
functional of <y?(x, t) and an odd functional of 7r(x, t). These observations allow us to deduce 
an ansatz for Qi that has the same structure as that in ()35|) : 

Qi = jjj dy dz M(xyz)7rx7ry7rz + JjJ ix dy dz A^x(yz)V5y7rx<^z, (68) 

where we have suppressed the time variable t in the fields and for brevity have indicated 
spatial dependences with subscripts. In ()68|1 the unknown functions M and have three 
arguments each. The function M is totally symmetric in its three arguments, and to em- 
phasize this symmetry we use the notation M(xyz); is symmetric under the interchange 
of the second and third arguments, and to emphasize this symmetry we write A^x(yz)- The 
functions M and A^ are like form factors because they describe the spatial distribution of 
the three-point interactions of the fields in Qi. We will see that the interaction of the fields 
is spatially nonlocal; this nonlocality is an intrinsic property of the operator C. 

We now proceed to determine M and A^. We substitute the ansatz (j68p into the commu- 
tator (j67|) and find after some algebra that two operator identities must hold: 



rfxfiyc/zrfw Ax(yz)G'^x¥'yV5w95z = -2 jdw(pl^, (69) 





c?x dy dz A'x(yz) (vTxVryV^z + fz'^^TTy) = 3 d:xidydzdw M^^y^^G^l^-Kyip^Ti^. (70) 
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By commuting ()69|) three times with vr, and ()70|) once with tt and twice with y?, we translate 
these two operator identities into two coupled partial differential equations for M and A^: 

{^? - V,2)iV,(y,) + (/x^ - Vy2)iVy(.,) + - V|)iV,(.y) = -65(x - y)5(x - z), (71) 

iVx{yz) + iVy(xz) = 3(/i2 - V,2)M(,y,). (72) 

To solve the system of coupled differential equations ()7H) and ()72|) . we Fourier transform 
to momentum space, denoting the D- dimensional Fourier transform of a function /x by 
/p = Jc?x /xC^P ''. Fourier transformation is effective here because it converts the differential 
equations (ffT|l and (ff^ into algebraic equations: 

iiVp(qr) + iiVq(pr) + iiVr(pq) = -6(27r)^(5(p + q + r) , (73) 
Gp Gq Gr 

iVp(qr) + iVq(pr) = iM(pqr), (74) 
Gr 

where Gp = (p^ + /i^)"^ 

Note that the right side of (j73|l contains the factor (5(p + q + r), which implies that 
the two three-point functions M and N conserve momentum. We thus introduce reduced 
representations of these vertex functions in which we have factored off the delta function: 

^(pqr) = (27r)^m(pqr)5(p + q + r), iVp(qr) = (27r)'^np(qr)5(p + q + r). (75) 

The functions in and h satisfy the following algebraic equations: 

i^p(qr) + i^q(pr) + i^r(pq) = "6, (76) 
Gp Gq Gr 

G'r^p(qr) + G'rraq(pr) = 3m(pqr). (77) 

There are two ways to solve these equations. A physically transparent but longer proce- 
dure making use of tree graphs is given in the Appendix ^ A shorter analytical approach 
is presented here. We begin by noting that the right side of ()77|) is totally symmetric in its 
indices. Thus, we can obtain two new equations by permuting the indices: 

GaP'Y>(cir) + G'qfir(pq) = 3m(pqr), G'pfiq(pr) + Gpnr(pq) = 3m(pqr). (78) 

We now have a sufficient number of algebraic equations, namely fj76l - ITH j) . to solve for in 
and n. The final results for M and A^ are 

M(pqr) = ^-^ - P - . ^(27r)^5(p + q + r), (79) 

P q ^p^r ^q^r — ^'-^'p'-^q'-^r 1,'-^ p T <-rq "T ^v) 
iVp(qr) = _ 6GpGqa(gpGr + Gpgq - Gq^r)^ _(2vr)^5(p + q + r) . (80) 

p w ; (^2(^2 A_ r<2 n2 i n2 n2 nn nn{n\n\n\ 
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As a check of these results we compare ()79|) and ()8U|) with ()36|) . which describes the case 
D = (quantum mechanics). When D = 0, we have just Gp = Substituting this 

expression for G into (f7^ and (jHUj) . we find that these equations reduce exactly to (IHUj). 

Next, we substitute Gp = (p^ + into and (jHUj) and use the inverse Fourier 

transform /x = (27r)^^ / (ip/pe"*^''' to express M and in coordinate space: 

fff dp dq dr 4e^*^-P-^yi-^^-'"(27r)^5(p + q + r) , . 

«i^y^> = JJJ . (81) 

fff dp d<i dr ^~"'"-'y"-" '{2i)''S{p + q + r) (q^ + - + (i^) 

= JJJ -Wr vC^) ■ 

where P(p, q, r) = p^ + q^ + r"^ - 2(p2q2 + q^r^ + r^p^) - 2^^{p'^ + q^ + r^) - S/i"^. We 
perform the r integral in (jHTjl and ()82|1 using the delta function and obtain 

4 ff gi(x-y)-p+«(x-z)-q 

^(xyz) = jjdvdq ^^p^^ , (83) 

where ©(p, q) = 4[p^q^ — (p ■ q)^] + 4/i^(p^ + P " Q + Q^) + S/i'* is positive, and 

iVx(yz) = 3 (Vy ■ V, + i/i^) M(xyz). (84) 



A. The (1 + l)-Dimensional Case 

For general D it is difficult to evaluate the double integral (jH^ in closed form. However, 
when D = 1 [(1 + l)-dimensional quantum field theory] we can evaluate the integral because 
the quartic terms in P(p, q) cancel. The evaluation procedure exploits the strict positivity 
of the denominator V{p, q), 

V{p, q) = 4:f/{p^ +pq + q^) + = 2ii'^[p^ + q^ + {p + qf] + 3/ > 0, 

to construct the one- dimensional integral identity T>^^ = dte^^^ {V > 0). This identity 
allows us to rewrite M(^^yz) as the triple integral 

M(,,,) = -^7^ [fdpdq rrfte-3'^'*e^(^-^)^'+*(^-^)^-^*(p'+^"?+^'). (85) 
l27rj JJ Jf-^Q 

Note that when D = 1 the variables x, ?/, z, p, g, and so on, are scalars and not vectors, so 
we no longer use boldface notation. 

To evaluate this integral we first complete the square in the q variable in the exponent 
and translate the q integration variable by g — g — | + ^^^j^- We then complete the square 
in the p integration variable and translate phy p ^ p — ^^^2^^ + ^^^et^^ ■ This gives 

= -j^^, jjdpdq j^^dte-'>^''-''^'-'''^'-P'l'^''^ (86) 

where p, which is totally symmetric in x, and z, is the positive square root of 

p' = \[{x-yf + {y-zf + {z-xf]. (87) 
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We now perform the scalings p — > p/v3t and q ^ q/ yM. The resuh is that the integral 
representing Mi^^yz) factors into three one-dimensional integrals: 



M, 



(xyz) 



-4/i"2(27r)-2l2-V2/2j, 



(88) 



where I is the Gaussian integral I = J dqe '^^ = a/tF and J = j^^Qdtt 



Finally, we use the integral representation [Ji 



/•oo 

Jt=o 



where Kq is the associated Bessel function. Thus, J = 2Ko(/ip). Combining the factors in 
we find that for a (1 + l)-dimensional quantum field theory (j83|l evaluates to 



M, 



1 



(xyz) 



Next we calculate using (jMjl . The result is 



Ko(/xp). 



(89) 



x{yz) 



3^3 
47r 



+- 



1^^ 



3(t/- 



Ko(/ip) + 



TT 



2p2 



KM 
pp 



5{x — y)S{x — z) 



(90) 



The mathematics underlying the solutions in (jH^ and ()9U|) is rather subtle and bears 
further discussion. First, it is important to mention that while we have expressed M and 

as functions of the three variables x, and translation invariance implies that these 
functions really depend on only two variables, say, the differences x — y and x — z. We 
therefore define the two variables rj and ( by 



X — ^{y + z) and 



In terms of these new variables we have 



(91) 



V3 



6{x — y)S{x — z) 



2 ^ 



2-Kp 



(92) 



The reason for introducing new variables and for emphasizing that we are working in the 
two-dimensional (r^, Q space is that in two-dimensional space the associated Bessel function 
^KQ^fip) is the Green's function: 



(/^' - Kc) 2^Ko(pp) = 6{v)6{0 = ^6{p). 



(93) 



This equation explains the appearance of the contact term (delta-function term) in the 
expression for in (jHOI)- The delta function is rotationally symmetric because the Green's 
function is rotationally symmetric and thus we can replace V^^^ in by d'^/dp^ + d/dp. 
Hence, we see that two derivatives of Ko(pp) give rise to a delta function: 



Kipp) = Ko(pp) - ^K'o(pp) - i<5(p). 



MP 



(94) 
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We have checked that M and in ()89|) and (jyup satisfy the partial differential equations 
()71|) and (f7^ . We can verify (f7^ by direct differentiation. To verify (f7H) we use the variables 
7] and C and take the indicated derivatives. The result is a combination of Ko(/ip), KQ(/ip), 
5'(p), and 5"(p) terms. The coefficients of all of these terms vanish except for the 
coefficient 5{p) term, and this coefficient reproduces exactly the right side of (ffT|l . 

Because our formulas for M and involve Bessel functions, we see clearly that Qi 
represents a nonlocal interaction of three fields. However, as associated Bessel functions 
decrease exponentially rapidly for large argument, the degree of nonlocality is small. 

VIII. QUANTUM FIELD THEORY WITH SEVERAL INTERACTING FIELDS 

The field theoretic calculations in Sec. lVIII can be extended to cubic quantum field theories 
having two and three interacting scalar fields. 

A. 'Pi<f'2 Theory 

We consider first the case of two scalar fields and y^x^'' whose dynamics is described 
by the Hamiltonian H = Hq"^^ + Hq^^ + eHi, where 

= \jd^ (7r(^))' + \ JJd^dy {G^Ji)-' (j = 1, 2) (95) 

and Hi = i jdn. (</3x '*)^V'x The Green's function Gxy is the solution to the equation 

(/^? - Vx^) = 5(x - y) (j = l,2). 

This quantum field theory is the analog of the quantum mechanical theory in (jllll . 
To determine C to order e we need to solve the operator equation 

= -2//i, (96) 

which is the two-field generalization of To find the solution to this equation we make 
an ansatz analogous to that in (|4ip for the operator Qi: 

where Mx(yz), A^yz, and ^x(yz) are unknown functions of three arguments each. As indicated 
by the parentheses, M and A^*^^-* are symmetric in their last two arguments. 
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We then substitute Qi into and use the following identities to perform the algebra: 



-1 



+2. jjjjd^d^dydz (GW)'^iV«,^f v^eVi^ 

+. jjljd^d^dydz {G^^iy' Niiy^v^v:' 

+1 [[[[dwd^dy dz {G<^1) M,(y,)7r«7r(i)(^t'^ • 



(97) 

Substituting (jHTj) into (j^Hjl gives an operator equation involving the unknown functions M, 
A^^^-*, and N^'^\ We then convert these operator equations into a system of three coupled 
partial differential equations by commuting with products of three fields: 



x(yz) 



25(x-y)5(x-z), 

.(2) 
x(yz)' 



(/i^ - V^) M.(y,) = iV«, + iV«y . (98) 

To solve this system of partial differential equations we perform Fourier transforms in the 
variables x, y, and z and obtain 



G 



(1) ^ Vqr ^ ^(1) ^ Vrq ^ ^(2) ' p(qr) 



-2(27r)^5(p + q + r) 



Grq 

Gp 



(99) 



Note that in (jHH) we have three linear algebraic equations in the four unknowns Mp(qr), 
Npqr, iVprq, and ^p(qr)- Thus, wc coustruct another equation from the second of the three 
equations in (j^Hl) by interchanging the momenta q and r: 



-(T)^p(qr)=iV«q + iVit)- 
(jrr 



(100) 



Solving the algebraic equations (jOHI) and (jlOOj) . we obtain the following expressions for 

M, iV(i), and iV^^); 



4 



p(qr) 



pqr 



p(qr) 



^^(p,q, r) 



(27r)^(5(p + q + r) 



2(p2 + q2-r2 + ;.2) 



I?(p,q, r) 
2[q^ + r^-p^ + 2^^-^i] 
P(p,q,r) 



(27r)^5(p + q + r), 

(27r)^5(p + q + r) 



;ioi) 
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where the denominator T>{p, q, r) is given by 

V{p, q, r) = + q^ + - 2 (p^q^ + ph^ + q^r^) + 2fil (p^ - q^ - r^) - ifijp^ + /i^ - i^ij^il 

Observe that if we set /ii = /i2 = 1 and take the quantum mechanical hmit D ^ 0, Eq. ()101|) 
reduces to (jl^ . 

Finally, we transform ()101|) back to coordinate space by calculating the inverse Fourier 
transforms. The integrals to be performed are triple D- dimensional integrals, but we can 
perform the integral over p by using the delta function. We get 

4 /"/" gj(x-y)-q+i(x-z)-r 

M.(y,) = -77^ lldqdr ^^7— ^ , (102) 



where P(q, r) = 4[q^r^ — (q ■ r)^] + 4/i^(q + r)^ — A^lq ■ r — /ig + 4:filfil, and 

Mill = [-V^ - Vy ■ V. + M.(y,), 



xyz 

x(yz) 



(103) 



We mention that for the special case D = 1 [quantum field theory in (1 + l)-dimensional 
space-time], the quartic terms in the denominator P(q, r) again drop out and it is possible to 
evaluate the integral in p02p in terms of Bessel functions using the the integration techniques 
described in Sec. IVllI For this special case we obtain 



x{yz) 



7r/i2A/4/i? - 



Ko(/^2p), 



(104) 



where 4/i|p^ = fil(2x — y — z)"^ + (4/if — fil){y — . The result in ()1()4|1 reduces to that in 



(jHn|l in the equal-mass case yUi = /i2 = Also, our results in D = 1 for A^iyl and A^^^^^-, are 



r(2) 



xyz 



47r/i2 
1 



Ko(/i2P) 



1 



^2 



2 _ .,2^{^- z){y- z) 
1 

2^il - {Afil - fil) 



(x- z){y- z) 



Ko(/i2p) 



5{x — y)S{x — z), 



^x(yz) 



V^pI-pI {x-z){y-z) 



47r/i2 
1 



Ko(/^2p) 



27rpnl^/4fil- fil _ 
1 



P2(4/i? -/ii) 



Alii - 2fii 



{x - z){y - z) 



P 



x{yz) 



KiP2p) 



5{x-y)5{x- z). (105) 



B. (pi(p2^3 Theory 



We now consider the case of three interacting scalar fields whose dynamics is described 
by the Hamiltonian 



(2) 



r(3) 



I5 



(106) 
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where Hq"^ is given in ()95p and Hi = i Jd:x. ip^'^ (p^^ (p^^ . This quantum field theory has the 
interesting property that its perturbative solution is finite for D < 3; there are no divergent 
graphs in less than three space-time dimensions. 

To find the operator C to leading order in e, we need to solve the operator equation 



r(2) 



(3) 



2) ^(1)^(2 
xyz z X 



+i 
+i 





+i 
+i 



(ix dy dz (^N^ 

dw dx dy dz (C^x 
dw c/x dy dz (G^x 
dxdydz (iVW^vrf 41 
dw dx dy dz (C^^ 
dw dx dy dz (C^], 





dxdydz (iV«,7r(i)7r(3 




+1 1 1 1 1 dw dx dy dz 
+i III dw dx dy dz 



-1 



From these equations we deduce the following system of differential equations: 

(/^? - Vx^)iVW, + (/i^ - V'^)N^t + (/^3 - V.^)iVi3)^ = -25(x - y)5(x - z) 





- ) Mxy. 


xyz ' yzx ' 


(/^? 




= A^(2) I 

yzx ' zxy ' 


(/^^ 




= Ar(3) + Ar(i) . 

zxy ' xyz 



(107) 



We introduce the ansatz 

Q, = jjjdxdydzN^ln^^^vf)^^)^ jjjdxdydzN^^y^^^f^^^) 

+ jjjdxdydzN^l7:^^~^^^)^^^^ + J^rfxrfy dzMx^.Tr^vrf vrf . (108) 
We then establish the following three results: 



(109) 



(110) 

The solutions for the unknown functions are as follows: Mxyz is given by the integral (j83p 
with the more general formula P(p, q) = 4[p^q^ — (p ■ q)^] + 4[/i^(q^ + P " q) + /^^(P^ + P ' 
q) — /i3P ■ q] + with 3m^ = 2/i^yU2 + '^lA^l^'i + '^lAl^'i " lA~ lA- "^^^ ^ coefficients 
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are expressed as derivatives acting on M: 



^xyz; 



xyz 5 
xyz- 



For the case D = 1 we have 



111) 



M, 



(xyz) 



Ko{mp), 



;il2) 



where 2m V = (/^? + 1^1- f^Di^ - y)^ + ifJ-l + IJ-l - l4)i.y + + I^Di^ -x?- We 
also have 



xyz 



\/3{x - y){x - z) 
47rp2 



Ko (mp) 



2npm 



yU2 + /^l-/^? V3(a;-?/)(a;-z) 



2p2 



-\/3 



5(x — — 2;), 



KM 

(113) 



(2) (3) 

and analogous expressions for Nxyz and Nxyi- 



IX. FINAL REMARKS 



We have introduced a new algebraic technique for constructing the operator C, which is 
required to define the positive-definite inner product of the Hilbert space in PT-symmetric 
quantum theories. Unlike the previously used analytical procedure for constructing C, which 
relies on the determination of all the energy eigenstates, the algebraic approach introduced 
here allows us to determine C directly from the operator form of the Hamiltonian. As a 
consequence, the approach extends naturally to quantum field theory. We have explicitly 
demonstrated the perturbative derivation of C both in quantum mechanics and quantum 
field theory for the case of cubic interactions. 

We hope to generalize the breakthrough reported in this paper to noncubic PT-symmetric 
quantum field theories, such as a —g(p^ theory. A —g(p'^ quantum field theory in four- 
dimensional space-time is a remarkable model because it has a positive spectrum, is renor- 
malizable, is asymptotically free [l^, and has a nonzero one-point Green's function Gi = (ip). 
Consequently, this theory may ultimately be useful in elucidating the dynamics of the Higgs 
sector of the standard model. 
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0^ 0^ 0^ 
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FIG. 1: Representation of -^p(qr) l|A4|) in terms of tree graphs. The solid lines represent factors 
of G and the dashed lines represent factors of G^^. The dots indicate points where momentum 
is flowing into the graph. The momentum in the graph is conserved and obeys the constraint 
p + q + r = 0. The numerical coefficients ai, 02, and multiply the graphs as shown in (|A4|) . 



APPENDIX A: GRAPHICAL SOLUTION TO EQUATIONS AND dTTj) 

Graphical methods can be used to solve the simultaneous equations (f7n|) and (f77|) for 
?7T.(pqr) and np(qr). We begin by defining the unknown function -Fp(qr) by 

?^p(qr) = Gp (-Fp{qr) " 2) . (Al) 

In terms of F, (f7n|) becomes 

Fpinr) + -P'q(pr) + -P'r(pq) = (A2) 

and (fTTjl becomes 

G'pG'q-Fp(qr) + G'qG'r-Fr(pq) " G'pG'r-Fp(qr) — G'qG'r-Fq(pr) = 2G'p (Gq — Gr) , (A3) 

where we have made use of fj74j) to eliminate rh. 

We can find an exact solution to ()A2|) and ()A3|) by using an iterative method. First, 
we construct a function that we call -f"p(qr)) which exactly solves ()A2|) and reproduces the 
structure on the right side of ()A3|) . An ansate that works is 

where the coefficients ai, 0:2, and are numerical parameters to be determined. 

The six terms in (IA4jl have a graphical representation (see Fig. in which we use a solid 
line to represent the Green's function Gp and a dashed line to represent its inverse G~^. 
The dots indicate momentum flowing into the graph. Note that the three momenta p, q, 
and r must satisfy a momentum conservation equation p + q + r = 0. 

The function F^^^ exactly solves ()A2|) if 

«! + 02 + "3 = 0, (A5) 
and F^^^ exactly reproduces the right side of ()A3|) if 



0^3 — 2a2 = 2. 



(A6) 
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FIG. 2: Representation of the ansatz -^p'(qr) ()A7|I in terms of tree graphs. The notation used in 
this figure is the same as that used in Fig. ^ The crosses in the graphs are vertices at which no 
momentum is flowing into or out of the graph. 

However, this choice for F^^^ also creates new terms on the right side of ()A3|1 that must be 
eliminated. To eliminate these additional terms we add to F^^^ a new ansatz, designated 
The form of F^^) is 

~(2) [Gl (gI Gl\ (gI G^\ Gl 

where the coefficients f3i, /?2; • • must be determined. The twelve terms in flA7|) are 
represented graphically in Fig. |21 

To determine the coefficients in F^'^'^ we require that the right side of ()A2|) vanish. This 
gives the three conditions 

A+/?2+/?3 = 0, 

+ = 0, 

Pe + 2/37 = 0. (A8) 

In addition, we require that the new terms arising on the right side of ()A3|) as a result of 
the ansatz F^^^ must be eliminated. This requirement gives the equations 

2a-i -ai- p5 = 0, 

ai - (32 + (33 - f37 = 0, 

a2 - (32 + (36 - (37 = 0, 

ai - a2 + (33 - (3e = 0. (A9) 
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FIG. 3: Representation of the ansatz -^p^qr) (|A10|) in terms of tree graphs. The notation used 
in this figure is the same as that used in Fig. |2 



One of these equations is redundant because if we subtract the third equation from the 
second, we obtain the fourth. 

We now continue the process. It is necessary to introduce a third ansatz in order to 
eliminate the new terms that appear on the right side of ()A3|) as a consequence of the ansatz 
F^'^\ This new ansatz is given by 



7(3) 
p(qr) 



, Gi\ Gl G%\ Gl Gl 



+74 



+77 




I(AIO) 



and it is represented graphically in Fig. |21 

To determine the nine numerical coefficients in F*^^) we require that the right side of ()A2|) 
vanish. This gives the three conditions 



71 + 72 + 73 
74 + 75 + 76 
77 + 78 + 79 



0, 



0. 



(All) 
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Also, we require that the new terms arising on the right side of ()A3|) as a result of the ansatz 
F*^^^ must be eliminated. This requirement gives the equations 





lb 


= 




- 76 


= 0, 


/?3 + /35-A 


-75 


= 0, 


pi -132 + 73 


-77 


= 0, 


A - 72 + 73 


-79 


= 0, 


/?2 - 72 + 77 


- 79 


= 0, 


/^T + 77 - 


-278 


= 0. 



(A12) 

Note that two of these equations are redundant; if we subtract the first from the second we 
obtain the third, and if we subtract the fourth from the fifth we obtain the sixth. 

It is clear that if we continue this process, we obtain more and more linear equations to 
solve. However, with each new ansatz the number of unknowns always exceeds the number 
of equations by one. For example, there are three a's but only two equations ()A5jl and ()A6jl : 
there are seven /3's but only six independent equations ()A8j) and ()A9jl : there are nine 7's but 
only eight independent equations ()A11|) and ()A12|) . Hence, with each new ansatz we have a 
kind of gauge freedom and we find it convenient to use this freedom to eliminate all graphs 
having split legs consisting of solid lines. That is, we choose Pe = = 0, 77 = 78 = 79 = 0, 
and so on. We can also choose ai = 0. Thus, for the next iteration our ansatz reads 

p(4) _ 
^p(qr) - 




Q2Q2 



(A13) 



The graphical representation of F^^^ is shown in Fig. |3] 

Requiring that the right side of ()A2|1 continue to vanish gives the three conditions 

5i + 52 + 53 = 0, 

Si + 65 + §6 = 0, 

2^7 + 58 = 0. (A14) 

Also, to eliminate the new terms arising on the right side of ()A3|) as a result of the ansatz 
must require that 



72 - 76 


-^4 + ^5 


= 0, 


73 + 75 


-76-54 


= 0, 


73 + 75 


-72-^5 


= 0, 


76 - 


- 274 + §7 


= 0, 


71 


-I2- h 


= 0, 




7i - "^i 


= 0, 


72 


- 5i + 53 


= 0. 



(A15) 
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FIG. 4: Representation of the ansatz -^p^qj.) iii l|A13|) in terms of tree graphs. The notation used 
in this figure is the same as that used in Fig. [21 Note that the number of hnear equations that 
determine the numerical coefficients is less than the number of coefficients. This allows us the 
freedom to exclude a priori all graphs having split legs consisting of solid lines. 

Note that two of these equations are redundant; if we subtract the first from the second we 
obtain the third ,and if we subtract the fifth from the sixth we obtain the seventh. Thus, in 
total there are eight 5's to be determined by eight independent equations. 

Let us now examine the solutions for the undetermined numerical coefficients. Solving 
the system of linear equations for the first five groups of these coefficients, we find that 





0, 


a2 = 
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"3 = 2, 
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53 = 2, 
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^5 = 14, 
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-56, 


57 


= 70, Ss 
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-2, 


£3 = 2, 


£4 = 


= 72, 


65 = 18, 


ee = 


-90, 


ej ■ 


= 252, eg 


eg 




-420. 





















-140, 
168, 
(A16) 

A brief examination of these coefficients shows that these numbers are all binomial coeffi- 
cients. In fact, by inspection we can now write down explicitly all terms in F^^\ F^'^\ F^^\ 
and so on, as the following double sum: 

" h G'i h ^2*/ Gl G^' + 1 - 2k} &■ 

It is easy to perform these sums, and the result for N is given in ()80p . It is remarkable that 
the solution for N is unique despite the arbitrary choice of "gauge" that we have made in 
solving the systems of algebraic equations. 
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We emphasize that all of the tree graphs that appear in this analysis have three external 
lines. This reflects the fact that the form factors M and that we are calculating represent 
the contribution of tree graphs to the three-point vertices. The next order in perturbation 
theory is proportional to e^, and the graphical representation of the result consists of all tree 
graphs having five external legs. 

[1] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 5243 (1998). 

[2] C. M. Bender and S. Boettcher and P. N. Meisinger, J. Math. Phys. 40, 2201 (1999). 

[3] E. Cahceti, S. Graffi, and M. Maioh, Comm. Math. Phys. 75, 51 (1980). 

[4] G. Levai and M. Znojil, J. Phys. A33, 7165 (2000). 

[5] B. Bagchi and C. Quesne, Phys. Lett. A300, 18 (2002). 

[6] Z. Ahmed, Phys. Lett. A294, 287 (2002); G. S. Japaridze, J. Phys. ASS, 1709 (2002); 
A. Mostafazadeh, J. Math. Phys. 43, 205 (2002); ibid. 43, 2814 (2002); D. T. Trinh, 
PhD Thesis, University of Nice-Sophia Antipolis (2002), and references therein. 

[7] P. Dorey, C. Dunning and R. Tateo, J. Phys. A 34 L391 (2001); ihid. 34, 5679 (2001). 

[8] A. Mostafazadeh, J. Math. Phys. 43, 3944 (2002). 

[9] C. M. Bender, D. C. Brody, and H. F. Jones, Phys. Rev. Lett. 89, 270402 (2002). 

[10] C. M. Bender, D. C. Brody, and H. F. Jones, Am. J. Phys. 71, 1095 (2003). 

[11] C. M. Bender, P. N. Meisinger, and Q. Wang, J. Phys. A 36, 1973 (2003). 

[12] C. M. Bender, G. V. Dunne, P. N. Meisinger, and M. Simsek, Phys. Lett. A 281, 
311-316 (2001). 

[13] C. M. Bender, J. Brod, A. T. Refig, and M. E. Renter, J. Phys. A: Math. Gen. (to 
be published). 

[14] R. Brower, M. Furman, and M. Moshe, Phys. Lett. B 76, 213 (1978). 
[15] M. E. Fisher, Phys. Rev. Lett. 40, 1610 (1978). 

[16] C. M. Bender, K. A. Milton, and V. M. Savage, Phys. Rev. D 62, 85001 (2000). 

[17] C. M. Bender and G. V. Dunne Phys. Rev. D 40, 2739 and 3504 (1989). 

[18] C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and 
Engineers, (McGraw-Hill, New York, 1978), Chap. 10. 

[19] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic, 
New York, 2000), p. 959. 



